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Drought is the most serious abiotic stress that hinders rice production under rainfed conditions. Breeding 
for deep rooting is a promising strategy to improve the root system architecture in shallow-rooting rice 
cultivars to avoid drought stress. We analysed the quantitative trait loci (QTLs) for the ratio of deep rooting 
(RDR) in three F 2 mapping populations derived from crosses between each of three shallow-rooting 
varieties ('ARC5955', 'Pinulupotl', and 'Tupa729') and a deep-rooting variety, 'Kinandang Patong'. In total, 
we detected five RDR QTLs on chromosomes 2, 4, and 6. In all three populations, QTLs on chromosome 4 
were found to be located at similar positions; they explained from 32.0% to 56.6% of the total RDR 
phenotypic variance. This suggests that one or more key genetic factors controlling the root growth angle in 
rice is located in this region of chromosome 4. 

Global climate change in recent years has caused serious drought damage to crop production, resulting in 
higher food prices, lower farm income, and more food insecurity. Therefore, strategies for the genetic 
enhancement of drought resistance in crop plants are essential for a stable and adequate food supply. Rice 
(Oryza sativa L.), a staple food for nearly half of the world's population, is susceptible to drought stress because of 
its shallow rooting compared with other cereal crops and therefore a limited capacity to extract water from deep 
soil layers 1,2 . Nevertheless, rice is frequently exposed to drought in 20% of its total growing area in Asia 3 . Many 
studies have suggested that a deep root system helps plants to avoid drought stress by extracting water from deep 
soil layers 4-6 . Typical upland rice has deeper rooting than lowland rice 7 . The deep root system in upland rice may 
contribute to its drought avoidance through enhanced water uptake 8 . Therefore, to ensure drought avoidance in 
shallow-rooting rice, introduction of genes responsible for deep rooting into such rice cultivars is considered as 
one of the most promising breeding strategies 6 . 

Deep rooting is a complex trait because the vertical root distribution is determined by a combination of the root 
growth angle and maximum root length in seminal and nodal roots of cereal crops 910 . Thus, the development of 
deep roots requires multiple genes that confer near- vertical root growth axis and increase root length along this 
axis. According to Courtois et al. 11 , a total of 103 quantitative trait loci (QTLs) for root length have been reported 
in rice. Recently, Obara et al. 12 fine-mapped qRL6.1, a QTL located on chromosome 6 and associated with root 
length in hydroponically grown rice seedlings. However, there have been few reports of mapped QTLs associated 
with root growth angle in rice. DROl, a QTL for deep rooting on chromosome 9 13 , and qSORl, a QTL for soil- 
surface rooting on chromosome 7 14 have been mapped. Both QTLs contribute considerably to the root growth 
angle. Cloning and characterization of DROl showed that it functions downstream of the auxin signalling 
pathway and controls the gravitropic curvature in rice roots 15 . Under upland conditions with drought stress, a 
near-isogenic line (NIL) with a functional allele of DROl introduced from the deep-rooting cultivar 'Kinandang 
Patong' has deeper roots and a significantly higher grain yield than the shallow- rooting parent variety 'IR64', 
which has a non-functional allele of DROl. Thus, DROl is the first gene associated with root system architecture 
(RSA) that has been shown to improve the ability to avoid drought. 

Sequence analysis of the DROl transcribed regions in 64 rice cultivars revealed six haplotypes 15 . Interestingly, 
some shallow- rooting cultivars grouped together with the cultivars carrying a functional allele of DROl such as 
'Kinandang Patong', suggesting the existence of undetected QTL(s) associated with root growth angle distinct 
from DROl . Availability of multiple genes that affect the root growth angle would be useful for breeding strategies 
using marker- assisted selection (MAS) to develop rice cultivars with high adaptability to drought. 

To identify new QTLs for root growth angle in rice, we selected three accessions (other than 'IR64') with the 
shallowest rooting from different DROl haplotype groups based on our previous study 15 . 'Tupa729' belongs to 
Group I together with 'Kinandang Patong'. 'ARC5955' and 'Pinulupotl' belong to Groups II and IV, respectively. 
We conducted QTL analyses for the root growth angle in three F 2 mapping populations derived from a cross 
between each shallow- rooting variety and 'Kinandang Patong' as a donor line with deep roots. 
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Figure 1 | Frequency distribution of the ratio of deep rooting in the three 
F 2 populations: AK-F2, 'ARC5955' X 'Kinandang Patong' (KP); PK-F2, 
'Pinulupotl' X 'Kinandang Patong'; TK-F2, 'Tupa729' x 'Kinandang 
Patong'. Vertical and horizontal lines above the bars indicate the mean and 
SD of each parental line. 

Results 

Phenotypic variations of the ratio of deep rooting (RDR) in the 
three mapping populations. Comparison of the RDR among the 
parental lines confirmed shallow rooting in 'ARC5955', 'Pinulupotl', 
and Tupa729' (ranging from 16.5 to 27.4%), whereas 'Kinandang 



Patong' showed deep rooting (approximately 80%) (Fig. 1). The 
RDR values in the three F 2 mapping populations from the crosses 
between 'Kinandang Patong' and 'ARC5955' (AK-F2), 'Pinulupotl' 
(PK-F2), and Tupa729' (TK-F2) were distributed between the values 
of the respective parental lines without considerable transgressive 
segregation (Fig. 1). The broad-sense heritability (h 2 B ) of RDR was 
79.9% in AK-F2, 73.6% in PK-F2, and 88.9% in TK-F2. 

Detection of QTLs for deeper rooting in AK-F2. Among 768 single- 
nucleotide polymorphism (SNP) markers, 346 (45.1%) were polymo- 
rphic between 'ARC5955' and 'Kinandang Patong'. To resolve the 
linkage gaps and to saturate the regions around the LOD support 
intervals for detected QTLs with DNA markers, we added 16 poly- 
morphic simple sequence repeat (SSR) markers in these regions. The 
AK-F2 linkage map, composed of 362 markers, covered almost the 
entire rice genome (Supplementary Fig. 1). The total map length was 
1533.5 cM, and the average distance between markers was 4.38 cM. 
Two QTLs for RDR on chromosomes 2 and 4 were detected based on 
an LOD threshold of 5.96 (Figs. 2 and 3a; Supplementary Fig. 1). 
The QTL on chromosome 4 had a large contribution to the 
phenotypic variance, explaining 32.0% of the total (Table 1). The 
mean RDRs of the lines homozygous for the 'Kinandang Patong' 
allele at the marker (RM6089) closest to the detected QTL were 
significantly higher than those of the lines homozygous for the 
'ARC5955' allele (Fig. 3b). The 'Kinandang Patong' allele at this 
QTL increased the RDR. The minor QTL on chromosome 2 ex- 
plained 9.5% of total phenotypic variance (Fig. 2; Table 1; Supple- 
mentary Fig. 1). A two-dimensional scan revealed no significant 
epistatic interaction in the whole genome in this population (Supple- 
mentary Fig. 2). 

Detection of QTLs for deeper rooting in PK-F2. Among 768 SNP 
markers, 184 (24.0%) were polymorphic between 'Tupa729' and 
'Kinandang Patong'. We also mapped 53 polymorphic SSR markers. 
The PK-F2 linkage map, composed of 237 markers, covered almost 
the entire rice genome (Supplementary Fig. 3). The total map length 
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Figure 2 | LOD score curves for the QTLs for the ratio of deep rooting in the three F 2 populations (see Fig. 1 for designations). Rectangles 
represent linkage maps with DNA marker positions shown as vertical lines. Chromosome numbers are indicated under each linkage map (short arms are 
on the left). Dotted lines indicate LOD thresholds (AK-F2, 5.96; PK-F2, 5.72; TK-F2, 5.46). 
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Figure 3 | Chromosomal positions and allelic effects of the QTLs for the ratio of deep rooting (RDR) on chromosome 4 detected in the three F 2 
populations (see Fig. 1 for designations), (a) The peaks of the LOD curves indicate the putative positions of RDR QTLs. Vertical lines in each linkage map 
indicate the genetic positions (cM) of DNA markers. Black bars on the linkage maps indicate 1 .8-LOD support intervals calculated by using the lodint 
function in R/qtP 6 . DNA markers are shown under the linkage maps; the numbers in parentheses indicate their physical map positions (Mb) on 
chromosome 4 of 'Nipponbare'. (b) Frequency distribution of RDR in three F 2 populations showing three genotype classes of the DNA markers closest to 
the RDR QTLs (RM6089 in AK-F2 and PK-F2; RM3916 in TK-F2). The inverted triangles indicate the mean, and the horizontal bars indicate SD of each 
allele. The same shading is used for triangles and corresponding bars. The means labelled with different letters differ significantly among the three alleles (P 
<0.05, Tukey's multiple comparison test). 
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was 1560.1 cM, and the average distance between markers was 
6.93 cM. Only one QTL (on chromosome 4) was found based on a 
LOD threshold of 5.72 (Figs. 2 and 3a; Supplementary Fig. 3). This 
QTL had a large contribution to the phenotypic variance, explaining 
36.2% of the total (Table 1). The mean RDRs of lines homozygous for 
the 'Kinandang Patong' allele at the marker (RM6089) closest to the 
detected QTL were significantly higher than those of the lines 
homozygous for the 'Pinulupotl' allele (Fig. 3b). The Kinandang 
Patong allele at this QTL increased the RDR. A two-dimensional 
scan revealed no significant epistatic interaction in the whole 
genome in this population (Supplementary Fig. 4). 

Detection of QTLs for deeper rooting in TK-F2. Among 768 SNP 
markers, 160 (20.8%) were polymorphic between 'Pinulupotl' and 
'Kinandang Patong'. We also mapped 64 polymorphic SSR markers. 
The TK-F2 linkage map, composed of 224 markers, covered almost 
the entire rice genome (Supplementary Fig. 5). The total map length 
was 1434.1 cM, and the average distance between markers was 
6.76 cM. Two QTLs were detected on chromosomes 4 and 6 based 
on an LOD threshold of 5.46 (Figs. 2 and 3a; Supplementary Fig. 5). 
The QTL on chromosome 4 had a very large contribution to the 
phenotypic variance, explaining 56.6% of the total (Table 1). The 
mean RDRs of lines homozygous for the 'Kinandang Patong' allele at 
the marker (RM3916) closest to the detected QTL were significantly 
higher than those of lines homozygous for the 'Tupa729' allele 
(Fig. 3b). The 'Kinandang Patong' allele at this QTL increased the 
RDR. A QTL detected on chromosome 6 showed a small genetic 
effect, explaining only 6.3% of total phenotypic variance (Fig. 2; 
Table 1; Supplementary Fig. 5). A two-dimensional scan revealed 
no significant epistatic interaction in the whole genome in this 
population (Supplementary Fig. 6). 

Discussion 

To identify new QTLs for root growth angle in rice, we selected three 
accessions with shallow rooting from different varietal groups of 
cultivated rice, and produced F 2 mapping populations derived from 
their crosses with 'Kinandang Patong'. We detected two QTLs for 
RDR in AK-F2, two in TK-F2, and one in PK-F2. These are novel 
QTLs for root growth angle because their chromosomal positions 
differed from those reported in our previous studies 13,14 . The three 
QTLs on chromosome 4, detected in the three populations, showed 
very high LOD values and a large contribution to the total phenotypic 
variance (Figs. 2 and 3; Table l).The physical positions of the mar- 
kers most closely linked to the LOD peaks of these QTLs indicated 
that they may correspond to the same locus. Therefore, we assume 
that one or more key genes responsible for root growth angle may be 
located in this chromosome region, and we designate this QTL as 
DEEPER ROOTING 2 (DR02). The 'Kinandang Patong' allele of 
DR02 increased the RDR in all three populations compared with 
the alleles from the three shallow-rooting parents. The 'Kinandang 
Patong' allele of DROl, detected in the 'IR64' X 'Kinandang Patong' 
recombinant inbred lines (IK-RILs), also resulted in deeper roots 
than the 'IR64' allele 13 . These results suggest that deep rooting of 
'Kinandang Patong' is largely controlled by DROl and DR02, 
although it is unknown whether these QTLs interact epistatically. 
To clarify the genetic interaction between them, further analysis 
using NILs for DROl and DR02 will be needed. 

We compared the relationship between the positions of DR02 and 
other root QTLs by using the QTL Annotation Rice Online Database 
(Q-TARO, http://qtaro.abr.affrc.go.jp/) 16 and the data summarized 
by Courtois et al. n . We found only one QTL for deep root ratio 
(percentage of the volume of deep roots [below 30 cm] relative to 
the total root volume), a parameter related to RDR, in the LOD 
support interval for DR02 (Fig. 3a; 27.74-30.69 Mb based on the 
'Nipponbare' genome sequence in the RAP database; http://rapdb. 
dna.affrc.go.jp/). This QTL was located between the SSR markers 



RM470 (28.28 Mb) and RM317 (29.25 Mb) in RILs derived from a 
cross between a lowland rice cultivar 'Zhenshan97' (indica) and an 
upland cultivar 'IRAT109' (tropical japonica)", but comparison of 
the marker positions alone is not sufficient to determine whether 
DR02 is related to this QTL. We also found that a QTL for root 
volume (qFSR4) ls had been fine-mapped near the candidate region 
of DR02. However, the candidate region of qFSR4 is located in the 
long arm at >31 Mb 18 , and thus is outside of the LOD support 
interval for DR02; this excludes qFSR4 as a candidate for DR02. 
According to the report of Courtois et al. n , many QTLs for other 
root traits (but not root growth angle) have also been detected in the 
putative region of DR02. Further study using a NIL for DR02 should 
be conducted to clarify whether DR02 has pleiotropic effects on 
other root morphology traits. For this purpose, we are currently 
developing a NIL for DR02 in each of the shallow-rooting varieties. 

We found no genes homologous to DROl in the DR02 LOD 
support interval. This suggests that DROl and DR02 are not 
structurally related. Many auxin-related genes such as CROWN 
ROOTLESS1 (CRL1)/ ADVENTITIOUS ROOTLESSV- 20 , OsPIDl 21 , 
CRL4/OsGNOMl 2X2i , and CRL5 24 are associated with root gravitrop- 
ism, which determines the root growth angle. Therefore, we checked 
the DR02 LOD support interval for the presence of auxin-related 
genes, but did not find any candidate genes in this region. We also 
compared the positions of DR02 and known root genes character- 
ized in mutants by using the Overview of functionally characterized 
Genes in Rice Online database (OGRO, http://qtaro.abr.affrc.go.jp/ 
ogro) 25 . We found three root-related genes, DWARF 176 26 , Oryza 
sativa indeterminate domain 10 27 , and GLR3.1 2 ", in the DR02 LOD 
support interval, but none of them has been reported to affect the 
root growth angle. Therefore, DR02 might be distinct from these 
known genes. 

Sequences of the DROl transcribed regions of 64 rice cultivars 
have been grouped into six haplotypes 15 . The three shallow-rooting 
varieties used in this study belong to different haplotype groups: 
Group I ('Tupa729' together with deep-rooting 'Kinandang 
Patong'), Group II ('ARC5955') and Group IV ('Pinulupotl'). The 
variety of Group I haplotype showed functional allele of DROl in the 
previous study 15 . Therefore, we expected to detect QTLs other than 
DROl in TK-F2. Indeed, we detected two new QTLs (but not DROl) 
accounting for 68.3% of the total RDR phenotypic variance in TK-F2 
(Table 1). This suggests that the differences in RDR between 
'Tupa729' and 'Kinandang Patong' are mostly explained by these 
two new QTLs. On the other hand, we expected that DROl might 
be detected in AK-F2 and PK-F2 because 'ARC5955' and 
'Pinulupotl' belong to different DROl haplotype groups than 
'Kinandang Patong'. However, DROl was not found in these map- 
ping populations. In primary mapping populations such as F 2 , 
segregation of major QTLs often masks the presence of minor 
QTLs and hampers their statistical detection 29 . We assume that 
DROl could not be detected in the present study because of a small 
difference in the genetic effects of the 'Kinandang Patong' DROl 
allele and those of 'ARC5955' or 'Pinulupotl'. This result also sug- 
gests that Group II and IV have functional allele of DROl. The 
detected QTLs in AK-F2 and PK-F2 contributed to only 39.8% and 
36.2% of the total RDR phenotypic variance, respectively. The h 2 B of 
RDR in IK-RILs was 77.7% 13 . In this study, AK-F2 and PK-F2 
showed similar h 2 B ; no significant epistatic interaction could be 
detected in either population. Therefore, we believe that other 
QTLs in addition to DR02 were associated with the RDR in these 
two mapping populations but were not detected in this analysis. For 
instance, we observed two additional LOD peaks in AK-F2 at lower 
LOD threshold (LOD = 3.40 between ah03001245 and ah03001697 
on chromosome 3; LOD = 3.78 at AD06014592 on chromosome 6) 
(Fig. 2; Supplementary Fig. 1). In PK-F2, we found four additional 
LOD peaks at lower LOD threshold (LOD = 4.80 at P0634_l and 
LOD = 4.09 at ah01003209 on chromosome 1; LOD = 4.24 between 
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AE02000044 and P0221_3 on chromosome 2; LOD = 5.08 between 
ad07001335 and ad07001571 on chromosome 7) (Fig. 2; Supple- 
mentary Fig. 3). Chromosome segment substitution lines are very 
powerful tools for detection of minor QTLs because they enable 
separation of minor QTLs from major QTLs located in other chro- 
mosome regions 29 . Such advanced materials will be useful to elucid- 
ate whether favorable alleles of DROl and other QTLs are present in 
AK-F2 and PK-F2. 

We conclude that DR02 is a promising QTL that can be used to 
enhance drought avoidance of shallow- rooting rice cultivars by con- 
trolling the RSA, because the 'Kinandang Patong' allele of DR02 
improves the RDR of cultivars in three different genetic back- 
grounds: aus, indica, and japonica. DR02 cloning will accelerate 
elucidation of its molecular function, and DR02 haplotype analysis 
using many rice accessions will enable breeding for a better RSA 
through MAS. Map-based cloning of DR02 is currently underway 
in our laboratory. 

Methods 

Plant materials. For QTL analysis, we developed three F 2 populations from crosses 
between 'Kinandang Patong 1 and (i) 'ARC5955' (138 plants; AK-F2), (ii) 'Pinulupotl' 
(134 plants; PK-F2), and (iii) Tupa729' (133 plants; TK-F2). 'ARC5955' (NIAS Acc. 
No. WRC35) is a landrace {axis) of Indian origin. 'Pinulupotl' (NIAS Acc. No. 
WRC24) is a landrace {indica) of Philippine origin. Tupa729' (NIAS Acc. No. 
WRC55) is a landrace (tropical japonica) of Bangladeshi origin. 'Kinandang Patong 5 
(International Rice Research Institute Acc. No. IRGC23364) is a landrace (tropical 
japonica) of Philippine origin. 

RDR measurement. RDR was evaluated quantitatively (as an index for root growth 
angle) based on a basket method that uses open stainless -steel mesh baskets (top 
diameter of 7.5 cm and depth of 5.0 cm; PROUD, Ushiku, Japan), as described 
previously 13 . We defined the RDR as the number of roots that penetrated the lower 
part of the mesh (i.e., >50 from the horizontal, centred on the stem of the rice plant) 
divided by the total number of roots that penetrated the whole mesh. Larger RDR 
values correspond to greater root growth angles. 

DNA marker analysis. We determined the genotypes of F 2 populations by using SNP 
markers selected from genome- wide SNP marker data 30,31 , and SSR markers selected 
on the basis of the data from the International Rice Genome Sequencing Project 32 . 
Total DNA was extracted from leaves by using the CTAB method 33 . The SNPs were 
detected from a set of 768 custom SNP panels by using a GoldenGate Genotyping 
Universal- 32 768-plex Assay Kit and the BeadStation 500 G system (both from 
Illumina, San Diego, CA, USA) according to the manufacturer's instructions. PCR 
amplification was performed in 5-ul reaction mixtures containing 0.5 ul (20 ng) 
DNA, 1.0 ul 5X PCR buffer, 0.1 ul 10 mM dNTPs, 0.025 ul (5 units) KAPA2G Fast 
DNA Polymerase (Kapa Biosystems, Boston, MA, USA), 0.125 ul of a mixture of 
forward and reverse primers (20 pM each), and 3.25 ul H 2 0. The PCR programme 
consisted of an initial denaturation for 1 min at 95°C; followed by 35 cycles of 10 s at 
95°C,10 sat55°C,andl s at 72°C; followed by a final extension for 30 sat72 D C.PCR 
products were separated by electrophoresis in 3% agarose gels (Type I agarose; Sigma- 
Aldrich, St. Louis, MO, USA) at 150 V for 120 min. 

Statistical and QTL analyses. h 2 B of RDR was estimated from the variances of the 
parental plants and the variance among the F 2 lines, and calculated as 

>>b[F2] = (Vp[V 2 ] - { Vp[Pi] + Vp[V 2 ]}/2)/Vp[V 2 ], 

where Vp[F 2 ] is the phenotypic variance among F 2 lines, and VpfPj and Vp[P 2 ] are 
the phenotypic variances of the parental plants. 

Linkage maps were constructed from the genotype data in MAPMAKER/EXP 3.0 
software 34 . Genetic distances were estimated by using the software's Kosambi map 
function 35 . QTL analysis was performed by using the composite interval mapping 
function of the R/qtl software 36 . The composite interval mapping threshold was based 
on the results of 1000 permutations at a 5% significance level 37 . The confidence 
intervals around each significant QTL peak were determined by 1.8-LOD support 
intervals 36 . The additive and dominant effects and the percentage of phenotypic 
variance explained by each QTL (R 2 ) at the maximum LOD score were estimated by 
using the sim.geno, makeqtl, and fitqtl functions in R/qtl 36 . Two-dimensional scans 
with a two-QTL model were conducted in R/qtl with the thresholds based on the 
results of 1000 permutations at a 5% significance level 36 . 
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